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Introduction 


This technical report was made as a supplementary part of the communication between the author 
and his supervisor, Professor Katsuhisa Furuta of the Tokyo Institute of Technology, during his 
time as a Ph.D. student there. 

Though the author call it a technical report, it is not always true that all of the material put in 
this report is new or novel. Since the author has not been very fussy in his choise of problems, the 
problems presented in here cover many thing and many of them are nothing more than interesting 
exercises. 


A system with absolute functions 


22" April 1998 
Figure 1 shows the phase plane of the system [Kha96] (Exercise 1.17(4)) 


L1 U+%2- 4X, (|v1| + |xa]), (CE) 
Gq = —2%, + %2— Xo (\x1| or |v2|) + U, 


when u = 0. 


Preliminary knowledge about the system 
The simultaneous equations 


Ly + L_ — £1 (|x 1| + |ro]) = 0, 
—2x1 + 2 — X2(|z1|+|z2|) =0 


are ill-formed, therefore cannot be solved for equilibrium points. 


Phase plot of the system 
Figure 1 and Figure 2 show the phase plane plot of Equation 1 
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Phase plane 


3 T T 


x2 
oO 
T 


Figure 1: Phase plane when there is no control. 
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Phase plot 


-1.5 -1 -0.5 0 0.5 1 1.5 


Figure 2: Phase plane when there is no control. toepy = 1.24sec simulated on Hayate. Initial 
conditions (x1, 2) = (4,2) and (i, —7), ¢ = —1.5, —1.4, —-1.3,...,1.5 
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Feedback of Signum function 


Figure 3 and Figure 4 show the phase plane when this system is subjected to an input 


x2 


u = sgn(x1 + 22). 


Phase plane 


Figure 3: 


Phase plane when the control input is u = sgn(x 1 + £2). 
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Phase plot 
2 T T T 


0.5- 


Figure 4: Phase plane when the control input is u = sgn(a1 + %2). topy = 1.49sec simulated on 
Hayate. Initial conditions (a1, 22) = (i, 7) and (i, —i), i = —1.5, —1.4,-1.8,...,1.5 
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Feedback of Signum function and states 


Figure 5 and Figure 6 Show the phase plane when 


u= —(@, + re) + sgn(a, + 22). (4) 


Phase plane 
3 T T 


x2 
Oo 
T 


Figure 5: Phase plane when the control is u = —(a, + 2) + sgn(a1 + 22). 


Discussion 


When there was no control. Figure 1 and Figure 2 shows that the equilibrium point at the origin is 
an unstable focus and, also that there is a limit cycle circling around the origin. The limit cycle is 
slightly larger in the x2 direction than in the x; direction. Every trajectory starting from an initial 
point other than the origin goes to and then stays on this limit cycle. 

When the control input is u = sgn(a#; + v2). Figure 3 and Figure 4 shows the discontinuity on 
the surface of the hyperplane. (s = 2; + 22 = 0) All trajectories still converge to the limitcycle, 
though the limit cycle is also distorted where it intersects the hyperplane. There is no node. Figure 
4 shows that every point which is on the hyperplane and inside the limit cycle acts as an unstable 
node. 

When the control input is u = —(x, + 2) + sgn(a1 + 22). Figure 5 and Figure 6 also shows the 
discontinuity on the hyperplane. Similar to Figure 4, Figure 6 shows that points on the hyperplane 
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Phase plot 


x OF 


Figure 6: Phase plane when the control is u = —(41 + 42) + sgn(a1 + 22). topy = 1.53sec simulated 
on Hayate. Initial conditions (x1, x2) = (i,7) and (i, —7), i = —1.5, —1.4, -1.3,...,1.5 


inside the limit cycle acts as an unstable node. In addition to Figure 4, Figure 6 has got two stable 
nodes. 
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Values of€ — nature of the origin 


—o <e<—2 stable node 
—2<e<0 stable focus 
e=0 centre 
Oe 2 unstable focus 
2<é€< 0 unstable node 
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Table 1: Characteristics of trajectories of a perturbed system with periodicity as depending on € 


Perturbed system with a periodic orbit 
Next consider the system |[Kha96] (Exercise 7.1(4)) 
t1 = Xo 
Eg = —2, +X, (1-2? - 23) tu, €>0. 
Some knowledge about the system 
By solving the simultaneous equations 


v2 = 


ats Sen a (1 =7} = 2) 


the equilibrium point is at (0,0). 
The Jacobian matrix at this equilibrium point is 


aft of 
Ox, 0x2 


0 1 
7 —l—2ery —2ey? + e(1 — 2? — y?) | 


x=0 


rl 


The Eigenvalues are then 5 + sv e*—4. Thus the nature of the equilibrium point at the origin 
depends upon the value of ¢. The characteristic of the trajectories around the equilibrium point is 


shown graphically in Figure 7 and as an intrepretation of Figure 7 in Table 1. 


Phase planes of various values of ¢ 


Draw the phase plane when u = 0 of this system as shown in Figure 8 (¢ = 1 x 107"), Figure 9 
(« = 1x 107°), Figure 10 (¢ = 0.001), Figure 11 (€ = 0.01), Figure 12 (« = 0.1), Figure 13 (« = 1) 


and Figure 14 (« = 5). The initial conditions used were either 
(41,22) € {5,-5} x {1 € {5 <1 < 5}},or 


or 
(x1, 22) = (i,é), (¢,-i), i= —1,—-0.8,-0.6,...,1. 
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Re1 and Re2 
Oo 
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| 
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Figure 7: Eigenvalues as a function of € of a perturbed system with periodicity 
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Phase plane 


4b 


-8 -6 -4 -2 0 2 4 6 


Figure 8: Phase plane. ¢ = 1 x 10-7 
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Phase plane 


4b 


Figure 9: Phase plane. ¢ = 1 x 107° 


Phase plane . 
: State plane (Epsilon = 0.001) 
T T 


Figure 10: Phase plane. ¢ = 0.001 No input. 
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Figure 11: Phase plane. ¢ = 0.01 No input. 


Phase plane 


State plane (Epsilon = 0.1) 
T T 


Figure 12: Phase plane. ¢ = 0.1 No input. 
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State plane (Epsilon = 1) 
1.5 T T T 


Figure 13: Phase plane. ¢ = 1 No input. 


‘ | | ; State plane (Epsilon = 5) 
al 
| ost 
*y | ost 
° | | oa 
: | | o.2b 
Teedltae | 
| ] g 0 
aac | -0.2 
yr | 0.4 
“SP | -0.6 
oi -0.8F 
4 2 
25 4 3 3 4 5 15 15 


Figure 14: Phase plane. ¢ = 5 No input. 
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Feedback of the signum function 


Introduce an input 
u = sgn (x1 + 22) 


and draw phase planes as shown in Figure 15 (¢ = 1 x 107"), Figure 16 (¢ = 0.001), Figure 17 
(¢ = 1), Figure 18 (€ = 5). 


Phase plane State plane (Epsilon = 10'-7)) 
T T T T 


Figure 15: Phase plane. ¢ = 1 x 10~" with VS input. 
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State plane (Epsilon = 0.001) 
T T T T T 


Figure 16: Phase plane. ¢ = 0.001 with VS input. 


Phase plane 


Figure 17: Phase plane. € = 1 with VS input. 
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Phase plane 


Figure 18: Phase plane. ¢ = 5 with VS input. 


Zh 
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With feedback of signum function and states variables 


When subjected to an input 
u= — (#1 + XZ) + sgn (x + 22) 


this system produces phase planes as shown in Figure 19 (¢ = 1 x 107’), Figure 20 (e = 0.001), 
Figure 21 (¢ = 1), Figure 23 (« = 5). 


State plane (Epsilon = 1E-7) 
T 


| Ze 


Figure 19: Phase plane. ¢ = 1 x 10~' with a VS input and feedback of state variables. 


An unstable case of positive state feedbacks plus a signum input 


Let the input be 
U=4%,+Xo+ sgn (x4 + £2) (8) 


which produce an unstable feedback system. Consider this case only briefly by showing the result 
in Figure 23 for various values of €. 


Discussion 


From Figure 8 to Figure 14 the trajectories spirals toward a limit cycle of radius one. But when ¢ 
becomes very small (¢ < 1) the rate which the trajectory moves toward this limit cycle becomes 
very small also, so that the result shown in Figure 8 and Figure 9 look like there are more than 
one limit cycle. Figure 12 to Figure 14 also shows that there are two hyper planes on either side 
of the limit cycle which the trajectory goes to and then slide toward the final destination. These 
hyperplanes become very distinct at ¢ = 5 as shown in Figure 14. 

When u = sgn (x1 + 22) the response depends on the value of ¢. With ¢ small xz; = +1 act as 
the two unstable foci as shown in Figure 15 (¢ = 107") and in Figure 16 (¢ = 0.001), while with 
larger € in the case of Figure 17 (« = 1) or Figure 18 (€ = 5) the state plane is divided into two 
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State plane (Epsilon = 0.001 
— T T 


| AS 


Figure 20: Phase plane. ¢ = 0.001 with a VS input and feedback of state variables. 


regions. Each region has got its own focus and the trajectory in each region converges to a small 
limit cycle which circle around the focus. 

When u = —21 — 42 + sgn (x1 + 22) Figure 19 to 21 show that there are two stable nodes at 
+0.5 for ¢ from very small upto « = 1. When ¢ = 5 (Figure 23) there is a limit cycle circling around 
both stable points which used to be nodes. Figure 21 and 23 shows hyperplanes which trajectories 
approach first and then slide along it toward the nodes (in the case of Figure 21) or toward the 
limit cycle (in the case of Figure 23.) 

Figure 23, though a stable feedback system, shows the effect of ¢ on the shape of the hyperplane 
and the trajectory which approaches it. Figure 23 from (a) to (f) are the result when ¢ becomes 
smaller and smaller. The trajectory in Figure 23 (f) looks like a sine wave. 
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State plane (Epsilon = 1 
T 


Figure 21: Phase plane. ¢ = 1 with a VS input and feedback of state variables. 


State plane (Epsilon = 5) 
1:5; T T T 


Figure 22: Phase plane. ¢ = 5 with a VS input and feedback of state variables. 
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Figure 23: Phase plane. (a) ¢ =5, (b) e =1, (c)e =0.1, (d) ¢ = 0.01, (e)e = 0.001, (f)e = 10~’, 
with a VS input and positive feedback of state variables. 
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A system with a periodic orbit 


21" April 1998 
Figure 24 shows the phase plane of the system [Kha96] (Exercise 7.1(1)) 


L1 = 9 (9) 
By = —a1 +2 (1— 3a} — 223) +4, (10) 


with u = 0. 


Some knowledge about this system 


The equilibrium point is the origin. (ie the solution of f(z) = 0 when « = f(x)) The Jacobian 
matrix of this system is 


230} 


emanate 
Ox 


0 1 Ot 
han eco eee al a) 


Then the Eigenvalues are 0.5 + 70.866 which means that the equilibrium point is an unstable focus. 


Phase plot of the system 
The phase plane of this system has got a periodic orbit. 


Phase plane 


Figure 24: Phase plane. No input. 


Technical report #1, presented to Professor K. Furuta by Kittisak Tiyapan 27 


With feedback of a signum function and state variables 


Figure 25 shows a phase plane of this system with a feedback of s.v. (state variables) and a VS 
control, ie. 


U=—2, — Lo + sgn(x1 + 22). 


Phase plane 
3 


Figure 25: Phase plane of the system with a feedback of state variables and a VS input. 


With feedback of a signum function 


Figure 26 shows a phase plane of this system with a VS input only, ie. 


u = sgn(x1 + 29). 


Discussion 


Figure 24 shows that the origin is an equilibrium point which acts as an unstable focus, and that 
there exists a limit cycle as well as two hyperplanes leading trajectories into this limit cycle. Both 
Figure 25 and Figure 26 show that there are two equilibrium points acting as nodes. It is interesting 


to note that in both of these two figures trajectories starting from the hyperplane will goto either 
one of these nodes. 


28 
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Figure 26: Phase plane of the system with a VS input only. 
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Another system with a periodic orbit 
Figure 27 shows a phase plane of the system [Kha96] (Exercise 7.1(2)) 


Ly = rH (12) 
Bo = =, bate — 2(ay + Qre)xe +-u, (13) 


with u = 0. 


Some knowledge about the system 


The solution of f(z) = 0 when « = f(x) gives the equilibrium point to be at the origin. The 
Jacobian matrix of this system can be found by 
0 1 


A= of 
Then the Eigenvalues are 0.5 + 70.866 which means that the equilibrium point is an unstable focus. 


0 1 
7 ae eae 


(Ox 


x=0 


Phase plot of the system 


Figure 27: Phase plane. No input. 
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Sample response with a variable structure input 


Figure 28 shows the phase plane when the system is subjected to 


u = sgn(x1 + 2). 


State plane 
T 


Figure 28: Phase plane with a VS input. 


Sample response with feed back of states and a variable structure input 


Figure 29 shows the phase plane when the system is subjected to 
u= —(@, + £2) + sgn(z, + £2) 


respectively. 
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Figure 29: Phase plane with a negative feedback of state variables and a VS input. 


ol 


a2 Tokyo Institute of Technology 


Figure 30 shows the phase plane of this system subjected to 
U=2,+ %2+sgn(r1 + 2X2). 


When u = x1 + %2 the phase plane is similar to Figure 30. 


Phase plane 


State plane 
T 


0.5- 


-0.5- 


Figure 30: Phase plane with a positive feedback of state variables and a VS input. 


Discussion 


Both Figure 27, 28 and 29 shows that there is a limit cycle as well as hyperplanes which lead to 
the limit cycle. Without an input Figure 27 shows that the equilibrium point at the origin is an 
unstable focus. Figure 28 shows that with a signum function input there are now three equilibrium 
points instead of one, the origin has changed into being an unstable node while the other two which 
are located at x; = £1,272 = 0 are now the unstable foci. 

Figure 28 also shows that the hyperplanes of the original system are still there while on the 
new hyperplane introduced via variable structure method (s = x, + 22 = 0) the trajectories are 
deflected. 

Figure 29 shows things similar to those which were described above for Figure 28. However for 
Figure 29 the deflection of trajectories at the hyperplane is more prominent and the two unstable 
foci are now at (+0.5, 0) instead of at (+1,0). (also the foci becomes more prominent in a way that 
trajectories wind around them more prominently) 

The unstable feedback shown in Figure 30 also has two hyperplanes. 
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DC motor 


20" April 1998 
An armature-controlled DC motor [Kha96] is shown here as Figure 31. The corresponding state 


re Ny in 

[| 
EO 
i 


Figure 31: Armature-controlled DC motor. 


equation is 


dw 

J ki (15) 
di 

5 ee eee «7 16 
a w—- Ritu (16) 


where 7, u, R, L are the armature current, voltage, resistance and inductance, J is the moment of 
inertia, w the angular speed. The constant excitation flux ¢ results in a torque ki and a back e.m_-f. 
kw 

Then the state equation of a DC motor with all the dimension made dimensionless is 


dw, 
1 = i, 17 
ht i (17) 
di,. ; 
Te = —W, — 1p + Ur, (18) 
where R TR r 
a) i u 
RE Gye Bat yh hk ech pa 


Tym is teh mechanical time constant, T, the electrical time constant and T;,, >> T.. Rewrite the state 
equation as 


E— = —Wp — ip + Up, (20) 
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where 
t ee Sn 


T, ©" 7, 7B 
Some knowledge about the DC motor system equation 
At the equilibrium point 


dy. di,. 
a =0, and d 


which gives the origin (7, = 0,w, =0) as the only equilibrium point. Moving « to the right-hand 
side in the second equation we have 


which gives the eigenvalues as 
1 


> (-1 ayia = 42)). 


Seeing that the eigenvalues depend upon the value of ¢, plot Figure 32 to show the values visually 
as € varies. The intrepretation of Figure 32 is Table 2. 


Re1 and Re2 
° 


Ll Ll Ll Ll 
-1 -0.8 -0.6 -0.4 -0.2 


0 1 1 1 1 
-1 -0.8 -0.6 0.4 -0.2 0 
Epsilon 


Figure 32: The eigenvalues as a function of ©. 


Response with various values of ¢ 


Figure 33 shows the phase plane of this system when ¢ = 5, Figure 34 « = 2, Figure 35 « = 1, 
Figure 36 ¢ = 0.1, Figure 37 ¢ = 0.01. 
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Values of ¢ Characteristic of the trajectories 


(—oco, 0) saddle point 

0 system degenerated 
(0, 0.25] stable node 
(0.25, oo) stable focus 


Table 2: Characteristic of the trajectories for DC motor system equation as function of € 


Phase plane of a DC Motor State plane (Epsilon = 5) 
T T T T T T 


Figure 33: Phase plane of a DC motor, ¢ = 5 


Phase plane of a DC Motor State plane (Epsilon = 2) 
T T T T T T 


Figure 34: Phase plane of a DC motor, ¢ = 2 
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Phase plane of a DC Motor State plane (Epsilon = 1) 
T T T T 4 T T T T T 


Figure 35: Phase plane of a DC motor, «= 1 


Phase plane of a DC Motor 
T T T 


State plane (Epsilon = 0.1) 
T T T 


Figure 36: Phase plane of a DC motor, ¢ =0.1 
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Phase plane of a DC Motor State plane (Epsilon = 0.01) 
T T T T T T 


Figure 37: Phase plane of a DC motor, ¢ = 0.01 
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Sample plots for various ¢ of the DC motor with states feedback and VS 

in the input 

Next with state feedbacks and a variable structure control, simulate the same system with an input 
Up = —lp — Wp + sen(i, + wr). 


The results are shown in Figure 38, Figure 39, Figure 40, Figure 41 and Figure 42. 


Phase plane of a DC Motor State plane (Epsilon = 5) 
T T T T T T 


Figure 38: Phase plane of a DC motor, ¢ = 5. With state feedbacks and a VS input. 


Discussion 
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Phase plane of a DC Motor State plane (Epsilon = 2) 
T T T T 4 T T T T T 


Figure 39: Phase plane of a DC motor, ¢ = 2. With state feedbacks and a VS input. 


Phase plane of a DC Motor State plane (Epsilon = 1) 
T T T T T T T 


Figure 40: Phase plane of a DC motor, ¢ = 1. With state feedbacks and a VS input. 
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Phase plane of a DC Motor State plane (Epsilon = 0.1) 
T T T T T T 


Figure 41: Phase plane of a DC motor, ¢ = 0.1. With state feedbacks and a VS input. 


Phase plane of a DC Motor State plane (Epsilon = 0.01) 
T T T T T T 


Figure 42: Phase plane of a DC motor, ¢ =0.01. With state feedbacks and a VS input. 


Technical report #1, presented to Professor kK. Furuta by Kittisak Tiyapan 41 


Van der Pol 


16" April 1998 
This time consider the Van der Pol system represented in a state equation form as 


Re SS (21) 


1 
es = -"2+2- 37 (22) 


Some knowledge about the Van der Pol system 


The only equilibrium point is at the origin. (z = 0, and —r+2- ae = 0) The Jacobian matrix 


1s 
0 1 
sl. a 3 
€ € € 


And the Jacobian matrix evaluated at the equilibrium point is 
0 i! 
OS Ay Wi csee | 
E € 


The eigenvalues of the system at the equilibrium point is then 
1 
— [1+ Vi —4¢| 
2é 


which can be shown graphically as Figure 43, the description of which is written in Table 3. 


| 
| 
Les 
m |e 
(Oe [ell seed 
a 


x=0,z=0 


a 
T 


° 


Re1 and Re2 
| 


L fl L 1 
1 -0.8 -0.6 -0.4 0.2 0 
Epsilon 


Figure 43: Exgenvalues movements when € varies for the Van der Pol system. 


Examples of response of Van der Pol with different values of ¢ 


Plot the graph of the state variables both in time domain and as a parametric plot for « = 1 (Figure 
44 and 45), « = 0.1 (Figure 46 and 47) and e = 1 (Figure 48 and 49). 
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Values of € nature of the origin 


(—oo, 0) saddle point 

0 degenerated system 
(0, 0.25] unstable node 
(0.25, oo) unstable focus 


Table 3: Characteristics of trajectories of the Van der Pol system as related to ©. 


i | | i it 
0 2 4 6 8 10 2 14 16 18 20 
Time in seconds 


-4 | Nl | | 


Figure 44: State variables of a Van der Pol system, ¢ = 1. 
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The phase behavior of the Van der Pol system 


Van der Pol (Epsilon = 1) 
T T T 


TO: 
, « Y)) : 


| 
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x 


(a) (b) 


25 -4 =3 


Figure 45: Phase portrait of a Van der Pol system when ¢ = 1. 


| l l | | 
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Figure 46: State variables of a Van der Pol system, ¢ = 0.1. 
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The phase behavior of the Van der Pol system 


Van der Pol (Epsilon = 0.1) 
T T T 


| l l | | 
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Figure 48: State variables of a Van der Pol system, ¢ = 0.001 
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The phase behavior of the Van der Pol system 


5 Van der Pol (Epsilon = 0.001) 
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Figure 49: Phase portrait of a Van der Pol system when ¢ = 0.001. 
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Van der Pol with feedback of states and variable structure in the input 


Next with a feed back 
T =—-—x —z+sgn(x + z) 


the phase plots of the Van der Pol system are shown in Figure 50, 51 and 52 fore = 1, ¢ = 0.1 and 
€ = 0.001 respectively. 


The phase behavior of the Van der Pol system 
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) oy, | 


= — 
7 -1.5- 
4 ol 
-2.5 


Figure 50: Phase portrait of a Van der Pol system with a variable structure control when ¢ = 1. 


Van der Pol (Epsilon = 1) 
T T T 


Discussion 


The Van der Pol equation has a phase response which still remains the same characteristic, having 
a limit cycle of more or less the a similar shape and appearance, Figure 45, 47, 49 are similar to 
Figure 50, 51, 52 respectively. The shape and appearance of the state plane trajectory depends on 
é. 
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The phase behavior of the Van der Pol system 


Van der Pol (Epsilon = 0.1) 
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2.5 


Figure 51: Phase portrait of a Van der Pol system with VS control when ¢ = 0.1. 


The phase behavior of the Van der Pol system 
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Figure 52: Phase portrait of a Van der Pol system when € = 0.001. 
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A second order system 


10" April 1998 
Consider the system 


g = Ar+Bu (23) 
y = Cr+Du (24) 
with 
—2 -l1 
ala a) 
1 
e-[o]: 
Gah | 
and D=0. 


Some knowledge about this system 


This system has got two poles at —1 +1 and one zero at —2. Figure 53 shows a root locus plot 
of this system. The only equilibrium point of the system is located at the origin because at the 
equilibrium point Ax = 0 which implies x; = 0 and x2 = 0. The corresponding eigenvalues of this 
system are 

—1+491.7321, 


and thus the origin is a stable focus for the trajectories. 
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Root locus 
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Figure 53: Root locus of the two dimensional system. 
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Sample response of the system with feedback of state variables and a 
signum function in the input 


Figure 54 shows state diagrams of this system with various types of state feedback together with a 
sign function. The inputs used to simulate Figure 54 are 


to 
io) 
aN 


L/) 


xX 0 XY 0 
-10 -2 \ 
-20 -4 ~~ 
-30 -6 
-4 -2 0 2 4 -5 5 
x1 x1 
20 6 


UA 


-4 -2 0 2 4 -4 -2 0 2 4 


Figure 54: From top down and from left to right, state diagrams of the system with a sign function 
together with a feedback of the two state variables, feed back of x,, feed back of x2 and feedback of 
neither one of the two state variables. 


u= Fxr+sen(CAzr + CBu) = Fx+sgn(—22+ wu), 


where Ff’ = [1 1], [1 0], [0 1] and {0 0] from top left first. In simulating to obtain Figure 54 algebraic 
loops and discontinuities were encountered. The simulation was stopped manually periodically to 
be able to obtain the result. 
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The state variables and the state plane of Equation 23 are shown in Figure 55. 


Figure 55: (a) The state variables (tsimutation = 10 sec, tcpu = 0.0100 sec, Initial conditions are 
(21,22) = (—7,—2)), and (b) the state plane (tsimutation = 10 sec, topu,average = 0.0077 sec, Initial 
conditions are (41,22) € {—10, 10} x {—10, 10}) 
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Let 
u = 3 sgn(x; + X2) 


be the input of the system and the result as shown in Figure 56 was obtained. 


Figure 56: (a) The state variables (tsimutation = 10 sec,tcpu = 0.0100 sec, Initial conditions are 
(241,03) = (—7,.—2)), and (b) the state plane (tsimiutation = 10sec; tepu average = 0:0102 sec, Initial 
conditions are (41,22) € {—10,10} x {—10, 10}), u = 3 sgn(x; + x2) 


Discussion 


Hyperplanes can be seen in Figure 54, though they might be somewhat strange in appearance since 
the hyperplane here is 
s(x) = —% +4, 


where uw is the input. Figure 55 shows that the system trajectory has a character of a stable focus 
with an equilibrium point at the origin. Figure 56 shows a division of the focus into two foci, system 
trajectory is then moving toward one of these two stable foci (which one it goes to depends on the 
initial starting point of the trajectory.) 
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Simple pendulum 


53 


Next consider a nonlinear problem of a simple pendulum. {[Kha96| (page 5) Let | be the length of 
the rod, m the mass of the bob and @ the angle subtended by the vertical axis through the pivot 
point and the rod . Assume that the rod is rigid and has zero mass. Let g be the gravitational 
constant and k be the coefficient on the bob moving through air the value of which is proportional 
to the speed of the bob. The equation of motion following Newton’s second law of motion is then 


mlé = —mgsin6 — klé 
The state equations of the system is 
L1 = 
t = ~9 sina, — —%2+ —, 
l m ml? 


providing x, = 6 and x9 = 0. 


Some knowledge about the system 
Let #, = #2 = 0 and solve the equations for the equilibrium points. We have 
to=0 “and -sin (a) = 0 
and therefore the equilibrium points are 
(ars Og ae EDs 
The Jacobian matrix is 
Of 0 1 
Ox | —¢cosz; —= 


and the Jacobian’s estimated at the equilibrum points are 


0 1 


x=(n7,0) 


which have the eigenvalues of 


1 


l 
| ta 22 Ay? 
ml i + 4/1 k? —4m igcos. 


2=nt 


As an example, let 


and A becomes 


, n=0,+1,42..... 


(25) 


(26) 
(27) 
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Eigenvalues of a simple pendulum 


T 
Nn 
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Figure 57: Eigenvalues of a simple pendulum (it’s real and imaginary components) 


as eigenvalues. The real and imaginary components of the eigenvalues are shown in Figure 57. 
From Equation 30 the term inside the square root is greater than zero when 


2 


l 
cos #1 < |—— 
: 4m?g 


= |K|, (31) 
here only the absolute value is necessary because in real systems 
Lk aU: 


but g can be either < 0 or > 0 depending on the convention used. According to Equation 31 there 
are two possibilities, that is 


1. Equation 31 is sometime true and sometime not. This is the case shown in Figure 57 where 
|| = 0.02. Another example is the case where 


1=21, k=167, ,g=—-9.8, m=0.7, (\K| =0.3). 


Figure 58 combines this case with Figure 57. 
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2. Equation 31 can never be true. This is the case demonstrated in Figure 59 where 


l=1, k=2, ,g=-98, m=0.14, (|K| =5.2). 


Eigenvalues of a simple pendulum 


Re1 and Re2 


Figure 58: Eigenvalues of a simple pendulum 


Depending on when and whether each of these two equations 


1 


l 
a ee yee Pe 
=| cyl k? —4m?lgcosx,;| < 0 (32) 


?k? —Am?lgcosz; > 0 (33) 


is true, the characteristic of the trajectory can be determined. The characteristic of the trajectory 
can then be summarized as shown in Table 4. 


Sample responses of a simple pendulum without an input 


Figure 60 show the state variables of this system. Here g = 9.8,! = 1, k = 1.5 and m = 3 were 
used in simulation. 


56 Tokyo Institute of Technology 


Eigenvalues of a simple pendulum 


Re1 and Re2 
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Figure 59: Exgenvalues of a simple pendulum 
Value of x1 Nature of the trajectory 
| (2nm — cos! im) ; (2nx + cos” sim unstable node 
x1 = 2nz + cos} iG there are equilibrium hyperplanes 
((2nm + cos! zx) P (2nm + cos! i)| stable node 
| (2m —cos-! i ; (2nx — cos! z{x;)) stable node 
((2(n —1)r+cos! jh) ,(2n7—) cos! 7) stable focus 


Table 4: Characteristic of trajectory for the simple pendulum 
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Figure 60: State variables of a simple pendulum with initial conditions x; = £2 = 8. No torque 1s 


applied. 
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Then with initial conditions (71,272) = (—10,10), (—6,10), (—2,10), (2,10), (6,10), (10,10), 
(—10, —10), (—6, —10), (—2, 10), (2,10), (6, -—10) and (10, —10) draw the state diagram of the 
system as shown in Figure 61. 


The phase behavior of the simple pendulum system 
30 T T 1 T 


-30 -20 -10 0 10 20 30 


Figure 61: State diagram of the simple pendulum with no input. 
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A simple pendulum with various kinds of input 


Figure 62 and Figure 63 are similar to Figure 60 and Figure 61 except that now there is an input 
(eae 


| 
0 5 10 15 20 25 30 
Time in seconds 


Figure 62: State variables of a simple pendulum with initial conditions x1 = x2 = 8. Constant input 
T =1. (Errata : the axis labels should be the other way around, ie the top figure should be labelled 
X_ while the bottom one x1) 
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The phase behavior of the simple pendulum system 
30 T T T T 
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Figure 63: State diagram of the simple pendulum with constant torque input. 
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Figure 64 and Figure 65 are also similar to Figure 60 and Figure 61 and now the input is a 
sinusoidal torque of a unit amplitude (ie T = sint). 


i 
0 5 10 15 20 25 30 
Time in seconds 


Figure 64: State variables of a simple pendulum with initial conditions x, = x2 = 8 and a sine 
input of a unit amplitude. (Errata : the axis labels should be the other way around, ie the top figure 
should be labelled x2 while the bottom one x1) 
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Figure 65: State diagram of the simple pendulum with a sinusoidal input. 
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Figure 66 and Figure 67 are also similar to Figure 60 and Figure 61 but with a feedback of state 
variables and a unit step input, ie 
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Figure 66: State variables of a simple pendulum with initial conditions x1 = x2 = 8 and a unit step 
input and a feedback of state variables. (Errata : the axis labels should be the other way around, ie 
the top figure should be labelled x2 while the bottom one x1) 
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Figure 67: State diagram of the simple pendulum with a unit step input and a feedback of state 
variables. 


Technical report #1, presented to Professor K. Furuta by Kittisak Tiyapan 65 


Figure 68 and Figure 69 shows the results when 
T = —2, — 2+ sgn(x2). 


The result shows a limit cycle. 
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Figure 68: State variables of a simple pendulum with initial conditions x, = x2 = 8 and a sign 
function input and a feedback of state variables. (Errata : the axis labels should be the other way 
around, ie the top figure should be labelled xz while the bottom one x1) 
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Figure 69: State diagram of the simple pendulum with a feedback of state variablesand a sign function 
input. 
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This time use a full state feed back together with a sign function of a hyperplane, that is 
T= 2+ 22+ sgn(xr1 + 22). 


The state variables as functions of time are shown in Figure 70 and the state diagram Figure 71. 
Figure 71 shows clearly sliding surfaces in the state plane. 
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Figure 70: State variables of a simple pendulum with initial conditions x; = ©2 = 8 and a hyperplane 
sign function input and a feedback of state variables. (Errata : the axis labels should be the other 
way around, ie the top figure should be labelled x2 while the bottom one x1) 


Discussion 


With no input applied, Figure 60 shows that the pendulum swings and the goes to an equilibrium. 
Figure 61 shows many equilibrium points all of which is the same position in space (though having 
different angles @) and having the angular velocity w = theta = 0. When a constant input is applied 
the pendulum turns around and around in circle varying the angular velocity periodically (Figure 
62 and Figure 63.) 

When the pendulum is subjected to a sinusoidal input it still goes to equilibrium as Figure 64 
and Figure 65 shows (though one might think that it should revolve on and on.) 

Feeding back state variables and one other function, Figure 67 whos the trajectory which goes 
to the equilibrium point at the origin, Figure 69 has got a limit cycle, while in Figure 71 there 
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The phase behavior of the simple pendulum system 


Figure 71: State diagram of the simple pendulum with a feedback of state variables and a hyperplane 
sign function input. 


are two equilibrium points representing a pendulum pushed and stay on either side of the bottom 
resting position. Figure 66, 68 and 70 are time-domain plots of the state variables for the cases of 
the state planes shown in Figure 67, 69 and 71 respectively. 
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A fourth order system 


24 April 1998 
To study the effect which the parameter ¢ has on the poles and zeros location and the shape of the 
root locus of a system with the following transfer function was considered. 


2(s + 5)(s? + 4s + 29) 


This system have got poles at —10, —1, —4+ 22 and —4 — 22; zeros at —5, —2 + 5i and —2 — 52; 
and the gain K = 2. The same system is then written in a state-space form. 


zt = Ax+Bu 
= Cxr+Du 
With 
—8 —)5 1 1.25 
4 QO 0 0 
a 0 0 —-11 -—2.5 |’ 
0 O 4 0 
0 
0 
as 2.82843 |’ 
0 
C= —2.82843 1.59099 0.70711 0.88388 
and 


D0: 


The root locus plot of this system is the plot of the root of the equation H(s) = 0 where 


H(s) = te and is shown in Figure 72. 


A singularly perturbed system 


Now consider the singular perturbation problem made by incorporating FE into the system equation 
above. 


Ex = Ax+Bu 
y = Cr+ Du 


Where EF = 


oo CO fF 
coor 
oro°; 
Nn OO ® 
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Figure 72: Root locus of the original problem. 


Root locus of the system with various values of ¢ 


Ther root locus plots of this system are shown in Figure 73, Figure 74, Figure 75 and Figure 76 for 
€ =.5,.1, .01 and .001 respectively. Figure 73 has got a transfer function 


2 (s + 10) (s? +45 +4 29) 


T(s) = 
(9) = G4 854 20) (5? + 115420)’ 
Figure 74 
T(s) = 25° + 108 s? + 458 s + 2900 
~ (s? + 85+ 20) (s? + 11s +100)’ 
Figure 75 
2s? + 1008 s? + 4058 s + 29000 
isl = 


(s? + 8s + 20) (s? + 11s + 1000) 


and Figure 76 
253 + 10008 s? + 40058 s + (2.9 x 10°) 


T(s) = 
(s) = 57+ 195? + 10108 5? + 802208 + (2 x 108) 
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Figure 73: Root locus when ¢ = .5. 


How « affects the system characteristic 


The value of ¢ which gives a pair of double poles on the real-axis was found by trials and errors 
to be « & 0.3306, at which point the pole location was approximately —5.5 (actual values are 
approximately —5.46 and —5.54) and the corresponding transfer function was 


2s? + 38.255? + 179s + 877.4 


T(s) = . 
(5) = 25 10s3 4 138.368? 4+ 4620 4 6081 


The other two complex poles remain at the same place all the time at —4 + 21. 

The gain margin remains all the time at oo while the phase margin remains all the time at 
deg 126.9141 at the frequency of 0.9730 rad/sec. The Bode plot, singular value plot, Nyquist plot 
and Nichols plot are shown in Figure 77, Figure 78, Figure 79 and Figure 80 respectively. 


As € approaches zero 


To observe the effect of ¢ on the system characteristic when ¢ — 0 show Figure 77, Figure 78, 
Figure 79 and Figure 80 again for ¢ very small. The Bode plot, singular value plot, Nyquist plot 
and Nichols plot are shown again in Figure 77, Figure 78, Figure 79 and Figure 80 respectively. 
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Figure 74: Root locus when « = .1. 
Some knowledge about this system 
4* April 1998 
Consider the system 
zg = Ar+Bu (34) 
y = Crt+Du (35) 
with 
—-§ -5 1 1.25 
4 0 0 0 
a 0 O —-11l —2.5 |’ 
Os: 20 4 0 
0 
0 
B= | 9.80843 | 
0 
C= —2.82843 1.59099 0.70711 0.88388 
and 


D0: 
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Figure 75: Root locus when ¢ = .01. 


From the system equation the only equlibrium point is the origin. (because Av = 0) The 
eigenvalues are 
—4+)2, -1 and —10. 


Response of the system with a feedback input 


Design a controller with an input signal 
u=Czrz+uU, 


where u; = sgn(CAx+CBu) and sgn(f) is called a sign function that returns the sign of the function 
f. Next simulate this system ten steps for six times, each with a different initial condition x(0). 
These (0) are (2 1 13 14)", [-2 —4 -5 —3]", [1 2 -3 —4]?,[-31 -2 5]?,(-132 —5]? 
and [2 —5 4 —3]?. Then draw state diagrams for some of the overall six parameter pairs. 

Figure 85 shows four such diagrams. Notice that there is a distinct sliding hyperplane in each 
diagram. 

With « = [5 5 5 5)]7,([-5 -—5 —5 —5]?, [55 —5 —5]7,[-5 —5 5 5]7,[-555 —5j? 
and [5 —5 —5 5]? as initial conditions do a similar simulation with the system with a 
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Figure 76: Root locus when ¢ = .001. 
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Bode Diagrams 
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Figure 77: Bode plot (¢ = 1 as a solid line, ¢ = 0.5 a dashed line, ¢ = 0.1 a dashdot line, both 
€ =0.01 and ¢ = 0.001 as dotted lines) 
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Figure 78: Singular value plot (¢ = 1 as a solid line, ¢ = 0.5 a dashed line, ¢ = 0.1 a dashdot line, 
both « = 0.01 and ¢ = 0.001 as dotted lines) 
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Nyquist Diagrams 
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Figure 79: Nyquist plot (« = 1 as a solid line, ¢ = 0.5 a dashed line, « = 0.1 a dashdot line, both 
€ = 0.01 and ¢ = 0.001 as dotted lines) 
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Figure 80: Nichols plot (¢ = 1 as a solid line, ¢ = 0.5 a dashed line, ¢ = 0.1 a dashdot line, « = 0.01 
circles and € = 0.001 cross marks) 
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Bode Diagrams 
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Figure 81: Bode plot (¢ = 1 as a solid line, ¢ = 0.001 a dashed line, € 


e€=107-' ande = 10~° as dotted lines) 


Frequency (rad/sec) 


10 


10-° a dashdot line, both 


80 Tokyo Institute of Technology 


Singular Values 
100 re 1 


-50 


Singular Values (dB) 


-100 


-150 


-200 | i i 
10° 10° 10° 
Frequency (rad/sec) 


Figure 82: Singular value plot (c = 1 as a solid line, ¢ = 0.001 a dashed line, ¢ = 10~° a dashdot 
line, both e = 10~" and ¢ = 107° as dotted lines) 
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Nyquist Diagrams 
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Figure 83: Nyquist plot (« = 1 as a solid line, ¢ = 0.001 a dashed line, ¢ = 10~° a dashdot line, 
€=10~' circles and ¢ = 10~° cross marks) 
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Figure 84: Nichols plot (¢ = 1 as a solid line, ¢ = 0.001 a dashed line, ¢ = 10~° a dashdot line, 
both « = 10-7 and e = 10~°® as dotted lines) 
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Figure 85: State diagrams of a system which « = 1. 
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Incorporation of ¢ 


The system equation 


¢ = E'Ar+E'Bu (36) 
= Cr+ Du. (37) 
Tt ih SOF 
Oe Te 
where F = 0010 
O20" -0.°e 


Some knowledge about the new system 


The equilibrium point still remains at the origin but the eigenvalues are affected by the ¢ introduced 


and are a a 
—44+)2, and — oa ae pve eave 


Sample plots of the new system with different ¢«’s 


Figure 86 shows a state diagram of a system which is equivalent to the system of Figure 85 but has 
a different set of initial conditions. 


Discussion 


The two poles which are not affected by the value of ¢ (the ones located at —4 + 21.) have the 
natural frequency w,, = 4.4721 and the damping factor 0.8944. The other two poles have the natural 
frequencies (w,,) and the damping factors (¢) as shown in Figure 89 and Figure 90 respectively. 

From Figure 72 to Figure 76 one can see that the two zeros at —2 + 51 remain at the same 
position what ever the value of ¢ may be while the other remain zero moves to the left away from 
the origin as ¢ + 0 at a high rate as shown in Figure 91. 

Figure 89 shows that w, — oo as ¢ — 0 whie Figure 90 shows that ¢ > 0 as e > 0. 

Both Figure 77 and Figure 78 show that the resonance freuency (w,) increases as € decreases. 
Figure 72 to Figure 76 show that the system is a minimum phase system Ve > 0, Figure 77, Figure 
78, Figure 81 and Figure 82 also show that this is true because the phase [Oga70| 


@ — deg —90(q — p) = deg —90(4 — 3) = deg—90, as woo, 


and the log-magnitude curve is 


—20(q — p) = —20(4 — 3) = —20 


asw —> OOo. 
decade 


Here p and q are the degrees of the numerator and denominator polynomials of the transfer function 
respectively. 
The Nyquist plot in Figure 79 shows that the system is always stable since [Oga70| 


Z=Nz=P=0 satisfies Z=N+P, 
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Figure 86: State diagrams of a system with e =1. 


where Z is the number of zeros of 1+G(s)H(s) in the right-half s plane, N the number of clockwise 
enciclements of the —1+ 70 point and P the number of poles of G(.S')H(s) in the right-half s plane. 
Both the Nichols plots of Figure 80 and Figure 84 show that at a given value of gain kK the peak 
or resonance value of M, 
Galea 
1+G 
increases as € decreases. 
With an input variable-structural in nature as shown from Figure 85 to Figure 88 the trajectories 
switch when coming to a switching hyperplane. 
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Figure 87: State diagrams of a system with e = 0.1. 
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Figure 88: State diagrams of a system with ¢ = 0.0001. 
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Figure 89: Natural frequency v.s. € 
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The Damping Factor 
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Figure 90: Damping factor v.s. € 
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Figure 91: The position of the moving zero v.8. € 
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